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We present a measurement of the azimuthal asymmetries of two charged pions in the inclusive 
process e + e _ —¥ ivkX, based on a data set of 62pb~ 4 at the center-of-mass energy of 3.65 GeV col¬ 
lected with the BESIII detector. These asymmetries can be attributed to the Collins fragmentation 
function. We observe a nonzero asymmetry, which increases with increasing pion momentum. As 
our energy scale is close to that of the existing semi-inclusive deep inelastic scattering experimental 
data, the measured asymmetries are important inputs for the global analysis of extracting the quark 
transversity distribution inside the nucleon and are valuable to explore the energy evolution of the 
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spin-dependent fragmentation function. 

PACS numbers: 13.88.+e, 13.66.Be, 13.87.Fh, 14.65.Bt 


The quark-hadron fragmentation process is 
parametrized with a fragmentation function (FF), 
which describes the probability that a hadron carry¬ 
ing a fraction of the parton energy is found in the 
hadronization debris of the fragmenting parton. The 
Collins FF, which considers the spin-dependent effects in 
fragmentation processes, was first discussed by Collins 
in Ref. lj. It connects the transverse quark spin with a 
measurable azimuthal asymmetry (the so-called Collins 
effect) in the distribution of hadronic fragments along 
the initial quark’s momentum. 

The measurement of the Collins FF provides an im¬ 
portant test in understanding strong interaction dynam¬ 
ics and thus is of fundamental interest in understanding 
QCD, the underlying theory of the strong interaction. 
Because of its chiral-odd nature, it needs to couple to an¬ 
other chiral-odd function, for instance, the transversity 
distribution M in semi-inclusive deep inelastic scatter¬ 
ing (SIDIS) or another Collins FF in e + e _ annihilations, 
to form accessible observables. The transversity distri¬ 
bution, which contributes to the nucleon transverse spin, 
corresponds to the tensor charge of the nucleon and is 
the least known leading-twist quark distribution func¬ 
tion. There have been several SIDIS measurements of 
this asymmetry from HERMES 0,0, COMPASS 0, and 
JLab 0 . Direct information on the Collins FF can be ob¬ 
tained from e + e _ annihilation experiments, as su gges ted 
in Ref. 0] . Measurements performed by the Belle [id IT3 | 
and BABAR 0 Collaborations give consistent nonzero 
asymmetries. Based on the universality of the involved 
functions in e + e _ and SIDIS [l4| experiments, global 
analyses min have been performed to simultaneous¬ 
ly extract the transversity and Collins FF. However, the 
e + e _ Collins asymmetries taken from Belle and BABAR 
correspond to considerably higher Q 2 (« 100 GeV 2 ) than 
the typical energy scale of the existing SIDIS data (most¬ 
ly 2—20GeV 2 ). Therefore, the energy evolution of the 
Collins FF at different Q 2 is a key factor to evaluate the 
transversity Recently, the treatment of the evolu¬ 

tion is developed in Refs [18H22I ]. which predict about a 
factor of 2 change in the observed asymmetries between 
BESIII energy and Belle and BABAR energy, but is not 
directly validated by experimental data. The BESIII ex¬ 
periment frit studies e + e _ annihilations at a moderate 
energy scale (4—20GeV 2 ). It is important to investigate 
the interesting feature of Collins FF at this energy scale, 
and the results can then be connected more directly to 
the SIDIS. Moreover, as emphasized in Ref. [l 8 (, with 
significantly lower Q 2 with respect to B factories, the re¬ 
sults will be crucial to explore the Q 2 evolution of the 
Collins FF and further the uncertainty of the extracted 
transversity, thus improving our understanding of both 
Collins FF and transversity. 

In this Letter, we present the measurement of az¬ 


imuthal asymmetries in hadron-hadron correlations for 
inclusive charged pion pair production e + e _ -A 7 T 7 tX, 
which can be attributed to the Collins effect. The anal¬ 
ysis is based on a data sample with an integrated lumi¬ 
nosity of 62pb _1 collected with the BESIII detector [23| 
at the center-of-mass energy = 3.65 GeV, where the 
energy is away from resonances. Compared to the ex¬ 
isting e + e~ data, in this measurement, only fragment¬ 
ing u, d, s quarks are involved. The results are free from 
charm contribution; as such, the combination with SIDIS 
data is more straightforward. The apparatus relevant 
to this work includes a main drift chamber (MDC), a 
time-of-flight system, and an electromagnetic calorime¬ 
ter (EMC). Details on the features and capabilities of 
the BESIII detector can be found in Refs. [23. l24j. 

Monte Carlo (MC) simulated events, which are pro¬ 
cessed with a full GEANT4-based [10 simulation of the 
BESIII detector, are used to optimize the event selection 
criteria and check for systematics. The MC samples for 
light quarks in e + e - —» qa (q = u, d , s ) processes are gen¬ 
erated by the LUARLW (20 package, which is based on the 
Lund model [HI, [28| . More MC samples including QED 
processes (e + e _ —> l + l~ (/ = e,/x,r), e + e _ —> 77 ), two 
photon fusion (e + e _ —> e + e~X), line-shape tail produc¬ 
tion of ip(?S) and, and the initial state radiative process 
e + e _ —> 7 J/ij) are analyzed to identify possible back¬ 
grounds. 

Taking into account the spin of the quark, the number 
density D'jJ for finding a spinless hadron h with trans¬ 
verse momentum produced from a transversely po¬ 
larized quark q with spin S q can be described in terms of 
the unpolarized FF D\ and the Collins FF H^ q at the 
leading twist 09j 
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where k denotes the direction of the initial quark q , 
z = 2 E h /Q denotes the fractional energy of the hadron 
relative to half of Q = \J s, and Mh is the hadron mass. 
The second term contains the Collins FF and depends on 
the spin orientation of the quark q , which leads to a sine 
modulation of the angle spanned by P^ and the plane 
normal to the quark spin. 

In hadron production in e + e _ —> qq events, the Collins 
effect can be observed when the fragments of the quark 
and antiquark are considered simultaneously. At y/s = 
3.65 GeV, due to the absent of the clear jet structure, 
there is no good way to estimate the q-q axis. However, 
the Collins asymmetries can be investigated with the az¬ 
imuthal angle (f >0 defined as the angle between the plane 
spanned by the beam axis and the momentum of the sec¬ 
ond hadron (P 2 ), and the plane spanned by the transverse 
momentum p t of the first hadron relative to the second 
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FIG. 1. The angle <j >o is defined as the angle between the plane 
spanned by the beam axis and the momentum of the second 
hadron (P 2 ) and the plane spanned by the transverse momen¬ 
tum pt of the first hadron relative to the second hadron. The 
angle 62 is the polar angle of the second hadron. 


hadron @ 113 , as shown in Fig. [Q 
The normalized dihadron yield is recorded as a func¬ 
tion of <j >0 and can be parametrized as a cos ([20o) + 6, 
with b referring to the term which is independent of 0o, 
and a can be written as @13 


a( 0 2 ,z 1 ,z 2 ) 


sin 2 0 2 T{Ht{z 1 )Ht{z 2 )/M 1 M 2 ) 

1 +cos 2 0 2 T [D 1 (zi)Di(z 2 )) 


( 2 ) 

where T denotes a convolution over the p t . The Mi and 
M 2 are the masses of the two hadrons, z\ and z 2 are 
their fractional energies, and d 2 is the polar angle of the 
second hadron with respect to the beam axis. D 1 and 
Hi denote FFs for antiquarks. 

We reconstruct charged tracks from hits in the MDC. 
We require the polar angle in the laboratory frame to sat¬ 
isfy |cos0| < 0.93, and the point of closest approach to 
the interaction vertex of e + e _ is required to be with¬ 
in 1 cm in the plane transverse to the beam line and 
within 10 cm along the beam axis. Particle identifi¬ 
cation for charged tracks is accomplished by combin¬ 
ing the measured energy loss (dU/dx) in the MDC and 
the flight time obtained from the time-of-flignt to de¬ 
termine a probability C(h = K, 7r, p , e) for each par¬ 
ticle ( h ) hypothesis. The tt ± (K ± ) candidates are re¬ 
quired to satisfy C(n)(C(K)) > 0.001, C(ir) > C(K) 
[C(K) > £(7r))and C(tt){C(K)) > C(p). Electrons are 
identified with the requirement C(e) > 0.001 and the 
ratio £(e)/(£(e) + £(n) + C(K)) > 0.8. Photons are 
reconstructed from isolated clusters in the EMC, whose 
energies are required to be larger than 25MeV in the 
EMC barrel region (|cos0| < 0.8) and 50MeV in end 
caps (0.84 < |cos0| < 0.92). It is required that the 
cluster timing delay from the reconstructed event start 
time does not exceed 700 ns in order to suppress elec¬ 
tronic noise and energy deposits unrelated to the event. 
To select inclusive e + e _ —> ttttX events, at least three 
charged tracks are required in order to strongly suppress 
two body decays. At least two of the charged tracks 
should be identified as pions. To suppress QED back¬ 
grounds with the final state t + t~ and unphysical back¬ 
grounds, e.g. beam-gas interactions, the visible energy 
in the detector, which is defined as the total energy of all 


reconstructed charged tracks and photons, is required to 
be larger than 1.5 GeV and no electron must be present 
in the event. Studies based on MC samples indicate that 
the backgrounds are suppressed to a negligible level, less 
than 2.5%. We select pion pairs with zi( 2 ) € [0.2, 0.9], 
where the lower bound is used to reduce pions originat¬ 
ed from resonance decays (mostly p, /), and the upper 
bound is used to reject two body decays. Compared to 
measurement at a higher energy scale 3: 3 , there is 
no clear jet event shape at BESIII which could help to 
separate the hadrons coming from different fragmenting 
(anti)quark. Instead, to select back-to-back pions, we 
require the opening angle of the two charged pion candi¬ 
dates to be larger than 120°. This requirement reduces 
the possibility that two pions come from the fragmenta¬ 
tion of the same quark. We label the two pions randomly 
as hi and /i 2 , and we use the momentum direction of h 2 
as reference axis. If more than two pions are present in an 
event, they are combined to each other, which means each 
pion is allowed to be assigned to different pion pairs. In 
the final event selection, 331 696 events survived, which 
provide 557 204 available charged pion pairs. 

We introduce the 20o normalized ratio, R = , 

where 7V(20o) is the dipion yield in each (20o) subdi¬ 
vision, and (No) is the averaged bin content. The nor¬ 
malized ratios are built for unlike sign r^), like sign 
(7r ± 7r ± ) and all pion pairs (7T7 t), defined as i? u , i? L and 
Rp , respectively, in which different combinations of fa¬ 
vored FFs and disfavored FFs are involved. A favored 
fragmentation process refers to the fragmentation of a 
quark into a hadron containing a valence quark of the 
same flavor, for example, u(d) —> 7r + , while the cor¬ 
responding u(d) —> tt~ is a disfavored process. Since 
the normalized ratio R is strongly affected by detector 
acceptance, we use double ratios RP/R Lt ' C ' 1 (UL and 
UC ratios) 0,13 to extract the azimuthal asymme¬ 
tries. The gluon radiation may induce a cos(20o) modu¬ 
lation according to Ref. 3, but it is highly suppressed 
at the BESIII energy scale and is independent of the 
charge of the pions. Through the double ratios, charge- 
independent instrumental effects cancel out, and QCD 
radiative effects are negligible at the first order, while 
the charge-dependent Collins asymmetries are kept. The 
double ratio follows the expression 


^L(cJ = A cos (20o) + B, (3) 

where A and B are free parameters. B should be consis¬ 
tent with unity, and A mainly contains the Collins effect. 
The Aul, Auc are used to denote the asymmetries for 
UL and UC ratios, respectively. 

The analysis is performed in bins of (z 1 , z 2 ), Pt and 
sin 2 0 2 /(l + cos 2 0 2 ). In (zi, z 2 ) bins, the boundaries are 
set at Zi= 0.2, 0.3, 0.5 and 0.9 (i = 1,2), where comple¬ 
mentary off-diagonal bins (z\, z 2 ) and (z 2 , z 1 ) are com¬ 
bined. In each bin, normalized rates i? u,L ’ c and double 
ratios R v /R h ’ C are evaluated. In Fig. [2J the distribu¬ 
tions of the double ratio / R L are shown as an exam- 
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FIG. 2. Double ratio R v /R h versus 2<% in the bin 21 € 
[0.3,0.5], 22 €[0.5, 0.9] (top) and bin 21 G [0.5,0.9], 22 € 
[0.5, 0.9] (bottom). The solid lines show the results of the fit. 



P t (GeV) 


FIG. 3. Asymmetries as a function of fractional energies 
( 21 , 22 ) (top) and pt (bottom) for the UL (dots) and UC 
(triangles) ratios, where the pt refers to the transverse mo¬ 
mentum of the first hadron relative to the second hadron, as 
shown in Fig. |T] In the top figure, the lower scales show the 
boundaries of the bins in 21 and 22 . Theoretical predictions 
from the authors of Ref. [l9j] are overlaid, where the hatched 
areas show the predicted bands. 


pie for two highest ( 21 , 22 ) bins with the fit results using 
Eq. ([3]). The asymmetry values (A) obtained from the 
fits are shown as a function of six symmetric ( 21 , 22 ) 
bins, p t and sin 2 0 2 /(l + cos 2 0 2 ) bins in Figs. [3] and [2 
respectively. The numerical results in each ( 21 , 22 ) and 
Pt bins are listed in Table H] 


< 0.1 

j- “I 

0.05 


0 


0.04 

7 7 

“ 0.02 
< 


0 



0.2 0.4 0.6 0.8 


sin 2 e 2 /(1+cos 2 0 2 ) 


FIG. 4. Asymmetries as a function of sin 2 ^ 2 /(l + cos 2 ^ 2 ) 
for UL (dots) and UC (triangles) ratios. Linear fits with the 
constant term being set to 0 (dashed line) or a free parameter 
(solid line) are shown. 


Several potential sources of systematic uncertainties 
are investigated and all systematic uncertainties are 
added in quadrature finally. An important test is the 
extraction of double ratios from MC samples, in which 
the Collins asymmetries are not included but radiative 
gluon and detector acceptance effects are taken into ac¬ 
count. In the MC samples, which is about 10 times of 
data statistics, double ratios are found to be consistent 
with 0 in all bins within statistical uncertainties. To test 
any potential smearing effects in the reconstruction pro¬ 
cess, MC samples are reweighted to produce generated 
asymmetries which vary in (0.02, 0.15) for UL ratios and 
(0.01, 0.08) for UC ratios in different bins. The recon¬ 
structed asymmetries are basically consistent with input, 
the differences between them, which range from 0.2% to 
48% for UL ratios and range from 2% to 31% for UC 
ratios relatively, are included in the systematic uncer¬ 
tainties. 

Additional possible contribution from gluon radiation 
can be examined in data by subtracting the normalized 
yields, f? u — R h ( c \ The subtraction method will cancel 
all the radiative terms, but the cancellation of the accep¬ 
tance effects may be incomplete. The differences between 
the asymmetries obtained with the subtraction method 
and the nominal results range from 0.001 to 0.01 for UL 
ratios and from 0.0 to 0.005 for UC ratios. These are 
assigned as absolute systematic uncertainties. 

The probability of misidentifying kaons as pions may 
introduce Kn pairs and KK pairs into the 7T7 t samples of 
interest. However, due to the much lower inclusive pro¬ 
duction cross section for charged kaons compared to pi¬ 
ons, the 7T7T asymmetry receives non-negligible contribu¬ 
tion only from the Ktt combination. We denote with A™ 
and A K * the corresponding Collins asymmetries in data. 
They can be obtained by unfolding the measurements of 
Ca. and <1, where A™ ea = (1 - f K ir)A™+f K TtA KlT , 
Amea. = ( x ~ Un)A Kn + f^A™, f Kn and f n7I are the 
MC-determined contamination fractions. Depending on 
the ( 21 , 22 ) bin, fn-* is found to range from 0.0% to 4.5% 
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TABLE I. Results of Aul and Auc in each ( 31 , 22 ) and pt bin. The uncertainties are statistical and systematic, respectively. 
The averages ( zt ), (pt) and are also given. 


Zl «"»• Z 2 

(zi) 

( 22 ) 

(Pt) (GeV) 

(sin 2 0 2 > 

Aul 

Auc 

(1+cos-^ 62 ) 

[0.2, 0.3] [0.2, 0.3] 

0.245 

0.245 

0.262 

0.589 

0.0128 ± 0.0085 ± 0.0114 

0.0050 ± 0.0038 ± 0.0017 

[0.2,0.3] [0.3, 0.5] 

0.311 

0.311 

0.329 

0.576 

0.0240 ± 0.0068 ± 0.0042 

0.0067 ± 0.0032 ± 0.0041 

[0.2,0.3] [0.5, 0.9] 

0.428 

0.426 

0.444 

0.572 

0.0281 ± 0.0131 ± 0.0077 

0.0136 ± 0.0064 ± 0.0029 

[0.3,0.5] [0.3, 0.5] 

0.379 

0.379 

0.388 

0.563 

0.0369 ± 0.0097 ± 0.0132 

0.0117 ± 0.0046 ± 0.0015 

[0.3,0.5] [0.5, 0.9] 

0.498 

0.499 

0.479 

0.564 

0.0518 ± 0.0120 ± 0.0049 

0.0217 ± 0.0056 ± 0.0046 

[0.5,0.9] [0.5, 0.9] 

0.625 

0.628 

0.499 

0.570 

0.1824 ± 0.0290 ± 0.0204 

0.0637 ± 0.0118 ± 0.0061 

Pt (GeV) 

(pt) (GeV) 

( 21 } 

( 22 ) 

(sin^ 0 2 ) 

Aul 

Auc 

(1 + COS^02> 

[0.00, 0.20] 

0.133 

0.291 

0.348 

0.574 

0.0122 ± 0.0093 ± 0.0021 

0.0044 ± 0.0043 ± 0.0006 

[0.20, 0.30] 

0.253 

0.285 

0.344 

0.579 

0.0279 ± 0.0081 ± 0.0034 

0.0100 ± 0.0038 ± 0.0016 

[0.30, 0.45] 

0.405 

0.327 

0.346 

0.570 

0.0241 ± 0.0072 ± 0.0025 

0.0090 ± 0.0031 ± 0.0026 

[0.45, 0.80] 

0.610 

0.453 

0.349 

0.571 

0.0516 ± 0.0087 ± 0.0040 

0.0211 ± 0.0049 ± 0.0019 

[0.80,1.40] 

0.923 

0.646 

0.334 

0.584 

0.0913 ± 0.0249 ± 0.0133 

0.0350 ± 0.0116 ± 0.0116 


and ranges from 0.1% to 35.4%, The errors on A K77 
are very large, and the changes in A 71 ’ 71 ’ from the nominal 
values are in (0.001, 0.005) for UL ratios and (0.0, 0.001) 
for UC ratios, and are assigned as systematic uncertainty. 

Additional higher harmonic terms (such as sin 20o and 
cos 40o) are also included in the fit function to validate 
the robustness of the fit. The changes of the value of the 
cosine asymmetries, which vary in (0.001, 0.009) for UL 
ratios and (0.0, 0.003) for UC ratios, are included in the 
systematic uncertainties. 

We have also verified null asymmetries for the double 
ratio of 7r + 7r + /7 t _ 7t _ pairs and for random combinations 
of pairs of tracks from different events. From these tests, 
no significant asymmetries are observed. The beam po¬ 
larization may contribute to the measured asymmetries. 
We study the angular distribution of the e + e _ p + p~ 
process, which is sensitive to beam polarization. No 
buildup of polarization is observed. 

Adding statistical and systematic uncertainties in 
quadrature, we observe significant, nonzero Collins asym¬ 
metries, as shown in Fig. [3] These asymmetries rise with 

I 


In summary, we perform a measurement of the 
azimuthal asymmetry in the inclusive production of 
charged pion pairs. Our results suggest nonzero asymme¬ 
try in the region of large fractional energy z , which can 
be attributed to the product of a quark and an antiquark 
Collins function. This is the first measurement of the 
Collins asymmetry at low energy scale ( Q 2 « 13GeV 2 ) 
in e + e _ annihilation. The observed asymmetry indicates 
a larger spin-dependent Collins effect than those at the 
higher energy scale from B factories ITol [bit . The results 
are of great importance to explore the Q 2 evolution of the 
Collins function and extract transversity distributions in 
nucleon. 


fractional energies and p t as expected theoretically [9] 
and seen in higher-energy e"*"e - experiments fici| - fl3l |. The 
predictions of authors of Ref. 0 , based on results from 
previous data and the energy evolution model, are also 
shown in Fig. [3j and are basically consistent with our re¬ 
sults. A direct comparison with higher-energy e + e _ da¬ 
ta is meaningless due to differing kinematics. However, 
asymmetries in our data are 1.5 times higher overall and 
higher by 0—2 sigma at points of comparable z and pt- 
The expected behavior of the Collins asymmetries as 
a function of sin 2 0 2 /(l + cos 2 0 2 ) is linear and vanishes 
at 62 =0, as formulated in Eq. 0- Thus, a linear fit 
is performed to the points in Fig. [I] with the constant 
term set to be 0 or left as a free parameter, which gives 
the reduced y 2 to be 2.3 or 2.8 for Aul and 1.7 or 1.9 
for Aye respectively. The significance for a zero offset is 
only about la for both Aul and Auc- 

The authors of the very recent paper Ref. 0 give the 
theoretical predictions for the BESIII energy scale, which 
are also shown in Fig. [3] Overall, our measured asym¬ 
metries are compatible with those predictions, except at 
the largest 2 interval. 
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